Strongly Correlated Superconductivity 



M. Capone\ M. Fabrizio^-^, C. Castellani\ E. Tosatti^'^ 
^ Dipartimento di Fisica, Universitd di Roma "La Sapienza" , 
and INFM Center for Statistical Mechanics and Complexity, Piazzale Aldo Moro 2, 1-00185 Roma, Italy 
^ International School for Advanced Studies (SISSA-ISAS), and INFM, 
Trieste-SISSA, Via Beirut 2-4, 1-34014, Trieste, Italy and 
^ International Center for Theoretical Physics (ICTP), P.O. Box 586, 1-34014 Trieste, Italy 

(Dated: February 1, 2008) 

High temperature superconductivity in doped Mott insulators such as the cuprates contradicts 
the conventional wisdom that electron repulsion is detrimental to superconductivity. Because doped 
fuUerene conductors are also strongly correlated, the recent discovery of high-critical-temperature, 
presumably s-wave, superconductivity in Ceo field efltect devices is even more puzzling. We examine 
a dynamical-mean-field solution of a model for electron doped fuUerenes which shows how strong 
correlations can indeed enhance superconductivity close to the Mott transition. We argue that the 
mechanism responsible for this enhancement could be common to a wider class of strongly correlated 
models, including those for cuprate superconductors. 



In conventional superconductors, the repulsive 
Coulomb interaction between electrons tends to oppose 
to phonon-mediated pairing, so that the actual critical 
temperature of the onset of superconductivity (Tc) 
decreases upon increasing electronic correlations. On 
the contrary, in the high-Tc superconductors, strong 
electron-electron correlations do not suppress super- 
conductivity, but rather seem to favor it because they 
are mostly poised on the brink of a repulsion-driven 
metal- Mott insulator transition (MIT). A recently 
developed approach capable of describing this transition 
is the so-called Dynamical Mean-Field Theory (DMFT), 
where spatial fluctuations are neglected, but the time- 
dependent quantum fluctuations are fully described [Q. 
As shown by DMFT, close to a Mott transition the 
large Coulomb repulsion U causes the effective metallic 
bandwidth W to be dramatically renormalized to a 
quasiparticle bandwidth = ZW , where Z ^1 \s the 
quasiparticle residue. This small effective bandwidth 
corresponds to an increased quasiparticle density of 
states at the Fermi level p = po/Z, which could at first 
sight be thought to enhance the attractive coupling 
A — pV and thus the critical temperature (po is the 
bare density of states at the Fermi energy per spin, V 
the pairing attraction). However, a decreasing Z does 
not automatically turn into an increase of A, because 
the pairing attraction V is itself renormalized down, by 
a factor Z^ within Migdal-Eliashberg theory, so that 
the increase of U finally depresses T^,, an effect further 
reinforced by a rising Coulomb pseudo-potential /i* . 

An enhancement instead of a decrease of with in- 
creasing Coulomb repulsion could take place if, on the 
contrary, the quasiparticle attraction did not get renor- 
malized by Z, and the repulsion U instead did. This 
is no doubt an appealing scenario, alas one which is at 
odds with all nai've expectations based on Landau Fermi- 
liquid theory. In this work we show that this scenario is 
actually viable. By solving a model for electron-doped 



Cgo , and closely examining the superconductivity arising 
there in proximity of the Mott transition, we find that 
correlations can indeed lead to a huge enhancement of 
phonon-driven superconductivity with respect to the un- 
correlated case. The analogy in the physics and even in 
the phase diagram of this fuUerene model as a function 
of decreasing bandwidth with that of cuprates for de- 
creasing doping draws a conceptual link between the two 
systems. 

The Hamiltonian describing this system is 
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where CR^i^ is the electron annihilation operator at site R 
in orbital i {i = 1, 2, 3) (the tiu level in Cgo is three-fold 
degenerate) with spin a, and nji — „nR,i„, where 

is the electron occupation number. 
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We also assumed for simplicity I^rr, = 5ijtRRi. We in- 
troduce the angular momentum density operators Li R = 

Ej,fc,^cJj,j<r^»jfcC^,fe<T' with i.jfc = -^e^jk proportional 
to the Levi-Civita tensor, and the spin density operators 
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the Pauli matrices. In terms of these operators, Hund's 
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term is Huund = 
liuR-if. 

The bare Jh is positive. However, in fuUerene, the 
Jahn- Teller coupling of electrons, and presumably also 
of holes ||], to the Hg molecular vibrations can reverse 
Hund's rules, favoring low spin and angular momen- 
tum ground states. We include this crucial electron- 
phonon effect by assuming Ju < 0, formally equiva- 
lent to treating the Jahn- Teller coupling in the antia- 
diabatic limit, where it can be shown to rcnormalize 
Jh Jh-SEjt/4: < 0, with Ejt the Jahn- Teller energy 
gain. The antiadiabatic approximation is justified for 
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fuUerene where vibron frequencies are as high as 0.1 eV, 
to be conipared with a correlation-narrowed quasiparti- 
cle bandwidth ZW, where the bare bandwidth 0.5 
eV and a quasiparticle residue Z <^\, due to a very large 
U/W. In any case, the neglect of retardation disfavors 
superconductivity, by preventing high-energy screening 
of the repulsion, hence overestimating fi^ . 

Following Ref. 3 we studied model (|l|) by DMFT(i1, 
varying U/W, at a fixed ratio Jh/U = —0.02 and integer 
filling (n) = 2 (or, equivalently, (n) =4). At weak cou- 
pling, U <^ W , model (|l|) describes a metal with three 
1/3-filled degenerate bands. If alone, the negative Jr 
would develop a superconducting s-wave order parameter 
Afl = Y^i=i cl? i|C^ jj^. With Jh = and considering ex- 
plicitly the electron-phonon coupling, Migdal-Eliashberg 
and DMFT calculations at relatively small U/W have 
well characterized this conventional, weakly or moder- 
ately correlated superconducting phase ||, ^, ^. How- 
ever, for our present Jh/U = —0.02 the effective super- 
conducting couphng A = 10po|'^ff|/3 — 0.2poU/3 (po is 
here the density of states per spin and band) is much 
smaller than the Coulomb pseudopotential /i* — poU, 
and weak coupling superconductivity is suppressed in fa- 
vor of a normal metal. At strong couphng, U » W, the 
system is a Mott insulator. Each site is occupied by two 
electrons which, since Jh < 0, form a spin and orbital 
singlet, as expected in a Mott- Jahn- Teller insulator |^). 
This state, a kind of on-site version of the Resonant Va- 
lence Bond (RVB) state 0, is nonmagnetic and has a gap 
to spin, orbital, and charge excitations. The transition 
between the metal and the strong coupling Mott insula- 
tor is however not direct, and a superconducting phase is 
known to intrude in between Q . The properties of this 
superconducting phase are, we now find, striking. 

Fig. 1 shows the superconducting gap A, obtained as 
the zero-frequency anomalous self-energy, compared with 
the hypothetical superconducting gap calculated in stan- 
dard Bardeen-Cooper-Schrieffer (BCS) theory at C/ = 0, 
and with the actual spin gap Agpin extracted as the edge 
of the main spectral peak in the dynamical spin suscep- 
tibility. 

Superconductivity is seen to arise suddenly out of the 
normal metal upon increasing repulsion above a critical 
value {U/W)c (here 0.8), and below the MIT (here at 
0.9). At {U/W)c the superconducting gap initially coin- 
cides with the spin gap, as in weak coupling BCS the- 
ory, but the two rapidly deviate. The larger spin gap 
merges with that of the Mott state; A instead reaches a 
peak value - a huge 1000 times the U = BCS gap cal- 
culated for the same pairing attraction - before falling 
again down to zero at the MIT. The large peak value 
of A is of the order of the maximum value which could 
be reached at J7 = if the bandwidth were comparable 
with \Jh\, which is also the condition to get the maxi- 
mum Tc for a fixed unretarded attraction |^ . In terms of 
a Landau Fermi-liquid description of the metallic phase, 
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FIG. 1: Superconducting gap A (squares) as a function of 
U/W = 50\Jh\/W. SC and MI stand for superconductor 
and Mott insulator, respectively. Also shown are the spin 
gap Aspin, reduced by a factor 20 (circles), and the BCS gap 
calculated as function of 50jj£r|/W at (7 = and multiplied 
by 1000 (long-dashed Hne). Gaps are normalized to X/po = 
10\ Jh\/3, which measures the pair attraction. 



this suggests that the quasiparticles close to the MIT 
have an effective bandwidth W* ^ ZW ^ \Jh\, and ex- 
perience an attraction of the very same order of mag- 
nitude. If this were indeed the case, the main effect of 
strong correlations would be to decrease W^ and thus in- 
crease the quasiparticle mass, leaving behind only a small 
residual quasiparticle repulsion. If the attractive vertex 
remained at the same time substantially unrcnormalized 
while Z — > 0, then the quasiparticle scattering amplitude 
would switch from repulsive at weak coupling to attrac- 
tive close to the MIT. A lack of renormalization of Jh 
is plausible, because Hund's coupling does not compete 
with U but rather benefits from it. In fact, U brings 
the system toward the atomic limit where Hund's rules 
are obeyed, whereas the metallic phase is where they are 
violated. There is here a similarity to the t — J model 
of cuprates, where J is also apparently unrcnormalized 
close to the insulator, as suggested by slave boson meth- 
ods jl^ and by numerical calculations ijlll . 

We found that this appealing, but thus far hypothetical 
scenario, is confirmed by a careful analysis of the metal- 
lic phase within Landau Fermi-liquid theory. DMFT en- 
ables a study of the normal metal even inside the super- 
conducting region, by preventing spontaneous breaking 
of the gauge symmetry in the self-consistency equations, 
and providing a full description of the quasiparticles and 
of their mutual interactions close to the Mott transition. 
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The Landau functional of the model, which possesses 
spin SU(2) and orbital 0(3) symmetry, contains here 
a multiplicity of Landau parameters f^^^\ gS{A) ^j^j 
/i^^-^^l). Defining i^-parameters F^(^) = 6po/^^^V^: 
G^(^) = 12po5^('^V^ and if^(^) = Apoh^^'^^/Z, di- 
mensionless quantities which measure the strength of the 
interactions between quasiparticles, the susceptibilities 
have the standard expression = ^J^' where x 
refers to the charge(spin) susceptibility for F — F^^'^\ 
and analogously for all the other orbital and spin-orbital 
(^QSiA) g^j^j ffS(A) parameters) susceptibilities. By cal- 
culating in DMFT the quasiparticle residue, Z, and all 
six susceptibilities, we obtain the F-parameters of the 
model as a function of U/W. 

Fig. 2 shows the decrease of Z in the metallic solu- 
tion on approaching the MIT. Superconductivity sets in 
at Z = Zcrit — 0.06, a very small value indeed. The 
charge compressibility decreases as a function of U/W, 
and vanishes at the MIT, consistent with the approach- 
ing incompressible insulator. The spin and all four or- 
bital and spin-orbital susceptibilities, which initially in- 
crease at small U/W (Stoner enhancement), turn around 
at U/W ~ 0.7, eventually vanishing at the MIT, consis- 
tent with a spin and orbital gap in the insulator. Ac- 
cordingly, monotonically increases from ^ U/W at 
weak coupling to infinity at the MIT, while the other 
parameters F^, G"^^"^-* and iJ^^^^) start off negative pro- 
portional to —U/W roughly until U ~ ZW, but then 
turns upward, cross zero, and finally diverge like at 
the MIT. 

This behavior of the Landau parameters in the metal- 
lic phase is at the root of the superconducting insta- 
bility, as is seen by calculating the quasiparticle pair 
s-wave scattering amplitude A. A has two contri- 
butions, Ai^i and Ai^j, describing singlet pair scat- 
tering from orbital i into the same or into another 
orbital, respectively. They are given by Ai^i ~ 
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and Ai^j = 
Fig. 2 shows 



At weak cou- 
= Jh, so that 
as the MIT is 



the pair amplitude A = Ai^.^ 
pling, A^^i = U + 4:Jh/^ > and A.^^j 
A — U + 10 Jh/5 is repulsive. However, 
approached and all F-parameters diverge, A —Z/2po, 
attractive and about equal to half of the quasiparticle 
bandwidth W*/2 — ZW/2, confirming our proposed sce- 
nario. In fact, the assumption of a quasiparticle repulsion 
renormalized by Z |p^ , and an unrenormalized attraction 
Jh would imply here A ~ ZU +10Jh /3. This simple ex- 
pression is seen to compare remarkably well with the true 
A up to {U/W)c leading to a very accurate estimate of 
0.067 for Zc„t- 

The crossover from weak to strong correlations occurs 
when the lower and upper Hubbard bands separate from 
each other, uncovering the quasiparticle resonance in the 




FIG. 2: Quasiparticle residue Z (triangles) as function of 
U/W. The vertical dashed line at U/W = 0.8 identifies the 
critical value above which the metallic solution may spon- 
taneously develop a superconducting order parameter. The 
quasiparticle scattering amplitude A (circles) and its heuris- 
tic approximation ZU + WJh /3 (long-dashed line) are also 
shown, both crossing zero at the metal-superconductor tran- 
sition. 



spectral function. This suggests a two-component de- 
scription of the model, similar to that used to analyze 
the MIT in terms of Kondo effect Q. We find that a 
two-component model describes very well the strongly 
correlated superconducting phase, in particular the prob- 
ability P{n) for a site to possess n electrons in the ground 
state. We calculated P(n) - the total weight of states 
with occupancy n in the ground state - for the super- 
conductor (Fig. 3 (B),(C),(D)), and that for the nearby 
Mott insulator Pins{n) (Fig. 3 (A)) to find that they are 
quite similar. In spite of an exceedingly large A, there is 
no evidence of preformed pairs or bipolarons in the super- 
conductor, as underlined by the strong steady peaking of 
P{n) around n = 2. 

By assuming the "two-component" form P{n) — 

ZP^p^\n) + (1 - Z)P^ns{n), where Z, P{n) and P^ns{n) 
are known, we extracted the quasiparticle probability 
distribution Pqp^'^ in the strongly correlated supercon- 
ductor between U/W = 0.8 and U/W = 0.9 (Fig. 3 
(A')(B')(C')(D')). It shows strong oscillations between 
even and odd n, as expected for a superconductor, with 
no major variations as a function of U/W even close 
to the Mott transition {U/W = 0.9), consistent with 
a weak to intermediate coupling superconductivity, im- 
plied by ^ ~ — W*. As a check, Pq^^^'^^\n) in the 
(metastable) non superconducting metal is also extracted 
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(Fig. 3 (A")(B")(C")) and found to be similar to that 
of a free Fermi liquid (Fig. 3 (D")), indicating almost 
free quasiparticles. The accuracy shown by this check is 
remarkable, as Pqp^^'^^'"^ is but a tiny fraction Z of 
P(n) and there is no free parameter. We conclude that, 
in the strongly correlated superconductor, free-fermion- 
like quasiparticles of weight Z become strongly paired 
while floating in a prevailing Mott insulator background. 
That background slows them down while taking away 
their Coulomb repulsion, but not their on-site (here Jahn- 
Teller originated) pair attraction. Somewhat similar to 
systems with spin-charge separation, the charge degrees 
of freedom are strongly renormalized close to the Mott 
transition, but the spin degrees of freedom - here includ- 
ing the pairing attraction - are not. The phase diagram 
of Fig. 1 for increasing U bears a remarkable similarity 
to that of cuprates for decreasing doping. We believe the 
superconductivity in the t — J model of cuprates to be in 
fact of a deeply similar origin - although the intersite an- 
tiferromagnetic interaction does of course introduce im- 
portant differences over our on-site Jh- The mechanism 
inducing singlet formation without competition with the 
Coulomb repulsion U is not far in spirit from Anderson's 
original RVB idea for cuprates Q. 

Coming to fuUerenes, we are only beginning to explore 
the full phase diagram and calculate Tc and other proper- 
ties for variable electron, and also hole doping. The possi- 
bility that superconductivity in these systems could be of 
the present, strongly correlated kind, seems real. Even if 
our solution is obtained for {n) = 2 or 4 (where supercon- 
ductivity has not yet been found), while the investigation 
of the (n) = 3 case (where superconductivity is actually 
observed) requires further work, we expect a very similar 
scenario also for the latter case. In particular, the chem- 
ically expanded (n) = 3 system (NH3)K3C6o is experi- 
mentally found to be insulating with low-spin (5 = 1/2), 
rather than high-spin (S* = 3/2) as expected for a regu- 
lar Mott state [Q. The low-spin is clearly of Jahn- Teller 
origin, indicating a Mott-Jahn- Teller insulator, exactly 
as in the (n) = 4 case of K4C6o- 

An important detail in fuUerenes is the actual value 
of the superconducting A. Strict electron-phonon cou- 
pling would yield a reahstic value of A ~ 0.8 — 
If alone, this large coupling would place fullerene super- 
conductors in the intermediate coupling regime already 
at C/ = 0. There, a further increase of A should not re- 
ally raise much Tc, or might even reduce it, in contrast 
with the well-known strong increase of Tc with increas- 
ing volume jist. An independent estimate of the effective 
pair attraction can be obtained by comparing the spin 
gap observed both in insulating K4C60 [|l6| and in su- 
perconducting K3C60 Aspin ~ 0.07 — 0.1 eV with 
that of our Jh < model, through Agpin — 5| Jjj-|[p^. 
We get in this way Jh ^ -0.02 eV ~ -0.02C/, hence 
A ~ 0.13, which is the tentative value adopted here. 
This large reduction of the effective A must, as explained 
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FIG. 3: Occupation probabilities for the particles, P{n), and 
for the quasiparticles in the superconducting and in the metal- 
lic phases, Pqp'~^\n) and Pqp'^'^^'"\n), respectively. 



above, be due to strong cancellation by the bare Hund's 
rule Jh 
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With a reduced A and large 
U , the only way to explain superconductivity and its in- 
creasing Tc with volume (increasing U/W), is to invoke 
the strongly correlated superconductivity we found just 
above {U /W)c- For further volume expansion, our model 
predicts as in Fig. 1 an eventual decline of Tc, and a Mott 
insulator for integer filling. Both features are observed 
in ammoniated compounds of the KaCeo family ]20[|. In 
the K4C60 family, conversely, U /W is above the MIT 
value, and we have Mott-Jahn- Teller insulators. Finally, 
we surmise that a similar strongly correlated supercon- 
ductivity, modified to account for the ci = 5 degener- 
acy of the hu hole states, should be relevant to the re- 
cently discovered C6o('^+) superconductivity pl| . Here 
a somewhat stronger electron-phonon coupling has been 
observed and calculated |Q, whereas the hole band- 
width and intra-molecular Coulomb repulsion are most 
likely similar to those of electrons. The enhanced su- 
perconductivity in CgoCHCls and CgoCHBra expanded 
lattices could result from the increase of U /W and of elec- 
tronic correlations, whereas the alternative explanation 
of a BCS-like increase in the density of states has been 
put in deep question by very recent results The full 
development of the theory of Ceo{n+) superconductivity 
is a task we reserve for the near future. 
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